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Abstract 

Property of charged fermion states is investigated in the quenched 
U(l) chiral Wilson- Yukawa model. Fitting the charged fermion propa- 
gator with a single hyperbolic cosine does not yield reliable results. On 
the other hand the behavior of the propagator including large lattice 
size dependence is well described by the large Wilson- Yukawa cou- 
pling expansion, providing strong evidence that no charged fermion 
state exists as an asymptotic particle in this model. 



1 Introduction 



The effort toward construction of chiral gauge theories on the lattice is pri- 
marily motivated by the fact that a full nonperturbative understanding of 
the standard SU(2) x U(l)y electroweak theory is required to answer a va- 
riety of unsolved questions such as the upper bound of the Higgs boson and 
fermion masses in the full theory and the magnitude of the baryon number 
violation in the standard model. The lattice approach provides a potentially 
useful tool for these purposes. However, there is a serious difficulty in this 
approach: because of fermion doubling |l|, each lattice fermion field yields 
16 fermion modes, half of which has one chirality and the other half the 
opposite chirality, so that the lattice theory is non-chiral. There are several 
attempts to avoid this difficulty. Among the proposals, the Wilson- Yukawa 
formulation p, H, ^ |[ has been extensively investigated during the last few 
years. In particular it has been found through numerical studies 0, |^, ^ 



and analytical work|IO|, |Il|, |T2|, [T^ that the fermion doublers decouple in the 
continuum limit if the Wilson- Yukawa coupling is sufficiently large. 

Recently it has been reported[0] in a quenched simulation of a U(l) 
chiral model that the charged fermion mass is approximately equal to the 
sum of the neutral fermion mass and the scalar boson mass on a 12^ x 16 
lattice. This relation suggests that the charged state is a scattering state of 
the neutral fermion and the scalar boson, and is not a bound state. Based 
on this interpretation it was claimed that charged fermion states do not exist 
in the symmetric phase of the global chiral model for strong Wilson- Yukawa 
coupling, and hence the Wilson- Yukawa formulation fails. 

It should be noted, however, that a large finite size effect is seen in the 
data of ref. for charged fermion propagators between 8^ x 16 and 12^ x 16 
lattices and that the relation among the masses does not hold on an 8^ x 16 
lattice. Thus much further analysis of the model, especially a systematic 
study of finite size effect, is needed to see if the conclusion is tenable. We 
have carried out such an analysis and we report on the results in this article. 

The organization of this article is as follows. In Sect. 2, we define the 
action of the U(l) chiral model we study and consider several bosonic and 
fermionic field operators, especially the two kinds of Dirac fermions, neutral 
and charged. We measure two-point correlation functions of these operators 
in the simulation whose parameters are also summarized. In Sect. 3 spectra in 
the broken phase are investigated. We extract masses of neutral and charged 
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fermions as well as those of scalar particles. In Sect. 4 we present and discuss 
the results of our Monte Carlo simulation in the symmetric phase in detail. 
Emphasis is placed on the fermionic charged state, which has been argued 



not to exist as a single particle in this phase 1 14]. Finite size effects for the 



charged fermion state are examined by means of the large Wilson- Yukawa 
coupling expansion. Our conclusions are summarized in Sect. 5. 



2 U(l) chiral model 
2.1 The action 

We investigate a chiral U(l) Higgs-fermion system in the Wilson- Yukawa 
formulation with the action given by 

Stotal ^boson ~l~ ^free ~l~ ^yukawa ~l~ ^wilson—yukawa 

— Sfjoson ~l~ S fermion^ (l) 

Sboson = -2Kj2Re{glgn+^), (2) 

Sfree = -^Y. (^nT/.V'n+M " i'n+^,'l^,'ipn) , (3) 



Syukawa 



Swilson—yukawa „ ^ ^ "^nPh (j]n+fi'^n+fi ~l~ 9n—^'4^n~fj, '^Qn^^'' 



_ 

-IjYl i^nPngn ii^n+f, + i^n-f, " '^i'n) ■ (5) 

Here gn is a non-linear U(l) scalar field satisfying gl^g^ = 1, and k, y and w are 
the hopping parameter of the boson, the Yukawa coupling constant, and the 
Wilson- Yukawa coupling constant respectively. The phase transition between 
the broken phase and the symmetric phase occurs at Kc = 0.149 in the 
quenched theory. The action Stotai is invariant under a global U(1)l x U(1)r 
transformation given by 

gn flLgn^R, ^L,R G U(1)l,R, 

V'„ ^ {QlPl + ^RPR)^n, 
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The scalar field g and the chiral fermion fields ipL and ipR are assigned the 
charges {Ql,Qr) = (1,-1), (1;0) and (0,1), respectively. 

We consider two kinds of Dirac fermion fields in our study defined by 

with the charge (0,1) and 

Cn={PL+9nPR)^n, = + oIPr) (7) 

with the charge (1,0). The fermion fields Nn and C„ are usually called 
"neutral" and "charged" in view of the eventual gauging of the global U(1)l 
symmetry. This distinction is not important since Nn will be the "charged" 
fermion if the U(1)r symmetry is gauged. The most important distinction 
between and C„ is that the Wilson- Yukawa coupling is reduced to an 
ordinary Wilson term and the Yukawa coupling to a mass term for the neutral 
fermion Nn- In fact one easily sees that 

3 fermion 7) ^ y ^n1 j^PhOn ^Qn+fi-^n+fi Qn—^-^n—fi) 

E NuI^iPr (A^n+M - Nn^^.) 
+yY,NnNn (8) 

n 

VO — 

-IT E (iVn+M + Nn^^. ' 2iV„) • 

Furthermore, the right-handed component of the neutral fermion has the 
shift symmetry[16]. Consequently, it is expected to appear as a massless 
particle at the point ?/ = 0. It has been shown ([0|-[|T3) that the physical 
fermion state of the Nn field appears in the continuum limit while doublers 
of the Nn are decoupled when the Wilson- Yukawa coupling w is large enough 
{w ~ 0(1) ). 

In terms of the charged fermion, on the other hand, the action is written 

as 

S fermion 'Z ^ ^ Cn'~i ^PhiCn+fi, Cn—fi) 
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+yY.CnCn (9) 
n 

X/ ^n9n (gli+nCn+ij, + Qn-nCn-^i — '^Qn^n) ■ 

Since w ~ 0(1) is needed for the doublers to decouple, the field Cn inter- 
acts strongly with the scalar field through the Wilson- Yukawa coupling. We 
investigate this strongly coupled system through the Monte Carlo simulation. 



2.2 Observables 

In order to extract the properties of the scalar state of the model, we consider 
three different correlation functions at zero spatial momentum defined by 

= X] < ^{Q,Q)^{n,t) >c, 
n 

D^if) = I] < 7^(0,0) m*) >c> (10) 

n 
n 

where (7„ = Re(g'„), 7r„ = Ini(gi„) and < >c means the connected part of 
correlation functions. 

In the symmetric phase (at k, < Kc = 0.149 in the quenched approxi- 
mation), only one massive scalar state is expected. Hence making a fit of 
form 

D^{t) = cosh[E^{t - T/2)] (11) 

with T the temporal lattice size should yield consistent values for the energy 
independent of the choice (!> = g,a,7T. 

The parameters and are related to the wave function renormal- 
ization constant and the mass parameter in the momentum-space- 
propagator through 

= 2A^sinh(E^)sinh(E^T/2), 
ml = 2cosh(E^)-2. (12) 
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In the broken phase (at k, < Kc) , the scalar field develops a vacuum expec- 
tation value and the symmetry U(1)l x U(1)r is broken down to U(1)l=r. For 
a finite volume, however, fiuctuation of the Nambu-Goldstone mode makes 
the expectation value of the scalar field vanish in absence of an external field. 
To avoid this problem, each scalar field configuration is rotated globally so 
that it satisfies J2n = 0. Consequently v =< g >=< a > and the mas- 
sive a propagator can be fitted by the same form as before. The form of 
the propagator for the Nambu-Goldstone boson is determined by the chiral 
perturbation theory [1T7| which gives 



^.(t) = A + 5((|-^r-3^), (13) 

where = AB/T. From the constraint J^t^Trit) = A and B are related 
by y4 = — i?/6T^. In our fit we take A and B as free parameters and check 
this relation for fitted results. 

For the neutral fermion, a general expression for the propagator follows 
from the shift symmetry ||16|| : 

Sf\p) = ij:i,^sm{p,){Ai{p)PL + A^{p)PR) + B{p), (14) 

A^ip) = (15) 
Bip) = y + wJ2{l-cosip^)). (16) 

The form of the function A^{p) can be evaluated with several analytical 



methods |TT|, |12|, |T^. The results are identical in the leading order, giving 

Aiip) = <gi^^gn>=z' = l/Z^. (17) 

In this approximation, the propagator at the zero spatial momentum takes 
the form 

^ < iV(o,o)%,t) >= PlIoSW) + PRloSm + S'J{t). 

n 

where 

^L,R(.. ^ ,L,R / cosh(E^(t - T/2)) cosh(ij;+(t - T/2)) \ 
^ ^ ' ^ cosh(E^T/2) cosh(E+T/2) )' 
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Mu^ _ AM ( smh(Ejv(t - T/2)) m + r(l - cosh E^y) 
V cosh(Ej^T/2) ' sinhE^ 



'S'at (^) — ^AT 

sinh(E^(t - T/2)) m + r(l-coshE 



cosh(E;^T/2) sinh E+ - N 



(18) 



with 



m = y 



VZ^Z^, r = wVZ^Z^, (19) 
" 2(m2 + 2rm + 1)1/2' (^0) 

^ 2(m2 + 2rm + 1)1/2' ^ ) 



/^+, m + r + V'HT'^ + 2rm + 1 
exp(^jv) = • (22) 

Here i?^ (-E'iv) is the mass of the physical fermion (the time doubler). With 
the replacement y — > y + 2w, the formula (|1^) applies to the (spatial) doubler 
mode at p = (vr, 0, 0). We denote its mass as End- 

We use the form (|TB|) taking and A^^''^^ as free parameters in order 
to fit Sj(f, S§ and Sff separately. The doubler propagator is fitted with the 
first term in (|T^ because its "time doubler" is much heavier. We will also 
compare data to the semi-analytical prediction (|18|) with the Monte Carlo 
estimation of as an input. 

For the charged fermion, the form of the propagator is not known. We 
found that the formula (|T^) does not fit the charged fermion propagator well. 



We therefore tried a single pole form||T^ for the fit given by 

S^'^'it) = A^'^ cosh(Ec(t - T/2)). (23) 

We extract the mass Ec using the data over a time interval tmin < ^ < 
T — tmin, varing tmin in order to see whether the mass is stable. 

2.3 Simulation 

Our simulation is performed in the quenched approximation on 8^ x 16 and 
12^ X 16 lattices. Scalar field configurations are generated by the 10 hit 
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Metropolis algorithm and the cluster algorithm of ref.[0 at k = 0.16 (broken 
phase) and at k = 0.145 (symmetric phase). See Table |l]for statistics. 

The fermion propagators are calculated by the conjugate gradient method 
with the point-source in the broken phase and the wall-source in the sym- 
metric phase. The anti-periodic (periodic) boundary condition in the time 
direction (the spatial directions ) is used. We fix the Wilson- Yukawa coupling 
at w = 1 and vary the Yukawa coupling y. The statistics depends on the 
simulation parameters. Typically, 2,400 (400) configurations separated by 
500 (200) sweeps are accumulated at k = 0.145(0.16) for the charged fermion 
propagator on an 8^ x 16 lattice. See Table ^ for further details. 

Errors are all estimated with the jack-knife method, using the bin size 
as listed in Table ^ and Table 0. For example, for bosonic observables on 
an 8^ X 16 lattice at k = 0.16, each bin consists of the average of successive 
20,000 (10,000) sweeps and the number of the bins is 4 (8) with the Metropolis 
(cluster) algorithm. 
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3 Results for the broken phase 



The TT propagator in the broken phase is shown in Fig.|I| together with the 
fit of form (0) expected from the chiral perturbation theory. Fit is excellent 
with A = -0.0218(52) and B = 23.7(2), which satisfies the constraint A = 
-5/6T2 = -0.0154(1) within 1.3a. 

The a propagator is shown in Fig.0 together with the fit of form (p^). 
The data points around t = T/2 lie above the fitting curve. We observed 
a similar behavior for the scalar field propagator obtained with the cluster 
algorithm |18|. This behavior may be related to the rotation procedure for 
calculating D^{t) in the broken phase. 

The masses and the wave function renormalization constants as well 
as the values of v and z"^ are listed in Table From the results we estimate 
the renormalized self coupling constant: 

= 3(^)X = 32.4(2). (24) 

av 

The tree-level unitarity bound for the 0(n) scalar model is Xr < 487r(n + 
1)-^ = 50.2(n = 2)[|1|]. 



In the broken phase we expect the relation Ej^r = Ec between the energy 
of the charged and the neutral fermion since they have the same quantum 
number under the residual U(1)l+r symmetry. This also can be seen from 
the following equation 



< CoCn >= V' < NoNn >+■■■. 

We fit the neutral fermion propagator with (|l^) and the charged fermion 
propagator with (p3|). Our results for En and Ec are plotted in Fig.^ by 



open circles and squares, (see Table § for numerical values.) They agree 
perfectly and also with the analytic prediction of the hopping parameter 



expansion [|T2|, |T3|, |T5[ drawn by dotted lines. Thus our results in the broken 
phase confirm the relation E^ = Eq assumed in refs.[^ (The relation 
En — Ec in the broken phase has been also observed in ref. using the 
momentum space propagators.) It is also seen in Fig.^ that the mass of the 
neutral fermion doubler End (solid circles) remains 0(1) in the y ^ limit. 
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4 Spectrum in the symmetric phase 



4.1 Scalar field 

The simulation has been performed with the parameters summarized in Table 
1^. The scalar propagator Dg{t) is shown in Fig.^ together with the fit. 
It is seen that finite-size effects are small. The result for the mass of the 
scalar field is summarized in Table ^ We have checked that these values 
are insensitive to the fitting range. Our value Eg = 0.415(1) on an 8^ x 16 
lattice obtained with the Metropolis algorithm is substantially larger than the 
valueQ Eg = 0.387(9) of ref. |]14|, which employed the cluster algorithm||18||. In 



order to make a direct comparison we generated 120,000 configurations with 
the same cluster algorithm and found Eg = 0.413(3) on an 8^ x 16 lattice, 
which perfectly agrees with our Metropolis value quoted above. Thus the 
discrepancy is not caused by a different choice of the algorithm. For those 
who may be interested in resolving this problem, we provide our data of 
Dg{t) in Table |. 



4.2 Neutral fermion 

The calculation of the neutral fermion propagator has been performed with 
the parameters in Table |^. The propagator data on an 8^ x 16 lattice are 
shown in Fig.|^. The propagators calculated through (|18]) (dotted curves in 
Fig.|^) agree with the data. Fitted values of the energy Ej^ are given in Table 
^ and are plotted as a function of y in Fig.|^ together with the prediction of the 
hopping parameter expansion. Again it is seen that the hopping parameter 
expansion agrees with the data. Finite size effects for E^ are very small. 



4.3 Charged fermion 

The propagator of the charged fermion S(j is shown in Fig.^. We observe a 
large finite size effect between 8^ x 16 and 12'^ x 16 lattices. As mentioned 
before we found that the propagators can not be fitted by the free fermion 
ansatz. We have therefore tried to extract the mass through a fit of the form 
( p3D over the range tmin < ^ < T — tmin- The value of Eq as a function 
of tmin is plotted in Fig.|. Since the value of Ec significantly changes with 

^This is calculated from the value nig = 0.390(9) quoted in ref.jT5| using (p^) 
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the variation of tmin, the charged fermion mass can not be extracted rehably 
by the fit. We are therefore not able to confirm the relation claimed in 
ref.|]T^ that Eq — + Eg on a 12^ x 16 lattice. We should comment that 



our propagator numerically differs from that of ref.[|14| by 3 to 5 standard 
deviations, exhibiting a faster decrease with t. (See Fig.|^.) Our data of Sq 
for both sizes are given in Table ^ The discrepancy probably originates from 
the difference of the scalar field configurations noted before. 

4.4 Interpretation of charged fermion state 

In order to understand the nature of the charged fermion state, we have to 
find a good explanation for both the t dependence of the rest energy Ec and 
the large finite size effect of the propagator. 



It has been reported in ref.[jT5[ that the large Wilson- Yukawa coupling 



expansion explains the data of ref.[14]. However, the agreement is only qual- 



itative. In addition our charged fermion propagator is different from that of 
ref.[n| and their results for the scalar mass and the wave function renormal- 



ization constant, which were used in the calculation of ref. ||T^ as an input, are 
also different from ours. Furthermore the large finite size effect of the charged 
fermion propagator was not discussed in terms of the large Wilson- Yukawa 



coupling expansion in ref. [0. We therefore think that a reanalysis of our 
data in terms of the large Wilson- Yukawa coupling expansion is neccesary. 

We have performed the expansion for the charged fermion propagator 
up to the next-to- leading order as in ref.[0. To this order the analytic 
expression for S^{t) is given by 

L L 



4 
1 

2M 



^(7)' I E sLv.it - ti)(v,,5L..(ti)) + wi'\t) I , (25) 



tl 



where a = 4/{Aw + y) is the expansion parameter and Vt/(t) = f{t + 1) 
f{t — l). The leading order contribution 

k 
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is the convolution of the neutral fermion propagator 



ikot sin(A;o) 

T ^ " {az/AY Ea sin2(A;A) + (1 - aw /A Ex cos(A;a))2 



ko 

and the scalar propagator 



The boundary conditions in the time direction imply /cqT = (2r;, + l)7r for the 
fermion field and PqT = 2n7r for the scalar field, where n = 0, 1, ■ ■ ■ , T — 1. 
In (p5D, W^\t) is the left-handed part of the 4-point- vertex contribution 
given by 

xG(^)(n,e,e + /i), 



with 



where G^^^ is the bosonic 4-point Green's function. Following ref. we take 



the mean- field result for G^^\ Similar formulae can be obtained for and 



Sq . See ref. for details. 



In order to compute the right handed side of (^Sj) we use the values of 
Zg, uig and determined in our numerical simulation results for and the 
simulation data for the two lattices are compared in Fig.p for y = 0.4 and 



in Fig.|TO|for y = 0.1. The agreement between the data and the prediction is 
excellent at y = 0.4 and is very good even at y = 0.1. 

The leading order contribution to Sq is the convolution of the neutral 
fermion propagator and the scalar propagator, each of which has very small 
lattice size dependence. It turns out, however, that their convolution has 
a large size dependence, which explains the size effect in our data. The 
next-to-leading order correction is found to be relatively small. The analytic 
predictions for S^i and S^^ also agree with the data very well. (See Fig.|Tl| 
and Fig.ll2D. 
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As a further check, we have performed the large Wilson- Yukawa coupling 
expansion of the charged fermion propagator in the broken phase. In this case 
the disconnected part has to be separated out from the two-point function 
of the scalar field, and the connected part is the sum of contributions from 
the massless Nambu-Goldstone boson and the massive real scalar particle. 
Therefore, we should add the following contribution 



k 



to (P3| ) and replace the scalar propagator Dg(t,p) with the sum of those for 
the Nambu-Goldstone boson and for the real scalar particle: 

1 ^r^' ™„+ 1 ^ jpijj Za- 



±y e^Pot + - y e 



T^^ 2Ea(1-cos(pa)) 2EA(l-cos(pA))+m2- 
Here means that po = is omitted in the sum for p = 0. The results up 



to the next-to-leading order are plotted in Fig.|T3|. The expansion works well 
also in the broken phase. 

Within the large Wilson- Yukawa coupling expansion it has been shown 
in ref. |TH] that the charged fermion is a scattering state of the scalar particle 



and the neutral fermion. The good agreement between the data and the 
result of the expansion therefore provides strong evidence that no charged 
fermion exists as a single particle state in the symmetric phase. 



5 Conclusion 

We have investigated the spectrum of fermions in the Wilson- Yukawa for- 
mulation through a quenched simulation of the U(l) chiral Wilson- Yukawa 
model. In the broken phase the neutral fermion state exists and the charged 
fermion operator generates the same state. In the symmetric phase the 
charged fermion mass can not be extracted by a single hyperbolic cosine 
fit even at large t and charged fermion propagators suffer a large finite size 
effect. We found that not only this large finite size effect but also the form of 
charged fermion propagators can be explained by the large Wilson- Yukawa 
coupling expansion. Within the expansion the charged fermion in the sym- 
metric phase is interpreted as a scattering state of the scalar boson and the 
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neutral fermion. We conclude that no charged fermion exists as a single par- 
ticle asymptotic state in the Wilson- Yukawa model in both phases. Only the 
scalar boson and the "neutral" Dirac fermion exist in this model. Although 



our data differ from those of ref. |TJ] and our reasoning which leads to our 
conclusion is quite different, we agree on the nature of the charged fermion 
state in the Wilson- Yukawa formulation. 

We can generalize this conclusion as follows. In more complicated chiral 
models with the Wilson- Yukawa coupling, we can always redefine fermionic 
field variables so that the Wilson- Yukawa coupling turns into a free Wilson 
mass term. Only those states generated by such field variables can appear 
as fermionic single particle states in the strong Wilson- Yukawa coupling re- 
gion. And these states are necessarily Dirac fermions with vector-like global 
charges. Their interaction with the scalar fields are weak in general. 
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Table captions 



Table 1. Simulation parameters for scalar fields. The symbol K represents 
1,000 configurations. 

Table 2. Simulation parameters for fermions. Metro, means Metropolis al- 
gorithm. Separation means the number of Monte Carlo sweeps between 
configurations. 

Table 3. Parameter and wave function renormalization constant for 
scalar fields and the value of v and z^. 

Table 4. Fermion masses in the broken phase on an 8^ x 16 lattice. The 
minimum value of t for fit is denoted by tmin- 

Table 5. Numerical value of the scalar field propagator Dg{t) in the sym- 
metric phase at K = 0.145 on an 8^ x 16 lattice. 

Table 6. Neutral fermion mass in the symmetric phase. 
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Figure captions 



Figure 1. a propagator D^{t) in the broken phase on an 8^ x 16 lattice. 
The solid line represents the fit. 

Figure 2. Same as Fig.l for the vr propagator DT^{t). 

Figure 3. Fermion masses as a function of y in the broken phase on an 
8^ X 16 lattice. Open circles (filled circles) represent the neutral fermion 
(neutral doubler) mass and open squares represent the charged fermion 
mass. The dotted line represents the prediction of the hopping param- 
eter expansion (HPE) for the neutral fermion mass. 

Figure 4. Propagator Dg{t) for the scalar field g in the symmetric phase on 
an 8^ X 16 lattice (filled circle) and a 12'^ x 16 lattice (open squares). 
The solid line represents the fit for the data on a 8^ x 16 lattice. 

Figure 5. Neutral fermion propagator S^{t) at y = 0.4 (filled circles) and 
y = 0.1 (open circles) in the symmetric phase on an 8^ x 16 lattice. Dot- 
ted lines represent the prediction of the hopping parameter expansion 
(HPE) for neutral fermion propagators. 

Figure 6. Neutral fermion masses as a function of y in the symmetric phase 
on an 8^ x 16 lattice (filled circles) and a 12^ x 16 lattice(open squares). 
Dotted line represents the prediction of the hopping parameter expan- 
sion (HPE) for the neutral fermion mass. 

Figure 7. Charged fermion propagator S(^{t) at = 0.4 in the symmetric 
phase on an 8^ x 16 lattice (open circles) and a 12^ x 16 lattice (open 
squares). Vertical bars connected by solid (dotted) lines represent the 
data of ref.|]n| for the same parameters on an 8^ x 16 (a 12^ x 16) 



lattice. 

Figure 8. Fitted value of Ec as a function of tmin in the symmetric phase 
aX y = 0.4 on an 8^ x 16 lattice (open circle) and 12^ x 16 lattice (filled 
circle). 

Figure 9. Charged fermion propagator S(^(t) at y = 0.4 in the symmetric 
phase on an 8^ x 16 lattice (filled circle) and 12^ x 16 lattice (open 



16 



circle), together with the prediction of the large w expansion including 
the next-to-leading order (solid lines). Dotted lines represent the con- 
volution of S^{t) and Dg{t), which is the leading order contribution in 
the expansion. 

Figure 10. Same as Fig. 9 at y = 0.1. 

Figure 11. Charged fermion propagators 5'^(t) at y = 0.4 in the symmetric 
phase on an 8'^ x 16 lattice (filled circle) and a 12^ x 16 lattice (open 
circle), together with the prediction by the large w expansion (solid 
hnes) . 

Figure 12. Same as Fig. 11 for S^{t). 

Figure 13. Charged fermion propagators, S(^{t), ai y = 0.3 in the broken 
phase. The meaning of symbols are the same as the Fig. 9. 
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Table 1: 





Broken 


Symmetric 


K 


0.16 


0.145 


size 


8=^ X 16 


8=^ X 16 


12^ X 16 


algorithm 


Metropolis 


Cluster 


Metropolis 


Cluster 


Metropolis 


T^sweeps 


80K 


80K 


1,200K 


240K 


900K 


#bin 


4 


8 


12 


24 


18 


bin size 


20K 


lOK 


lOOK 


lOK 


50K 



Table 2: 





Broken 


Symmetric 




0.16 (Metro.) 


0.145(Metro.) 


size 


8^ X 16 


8^ X 16 


12^ X 16 


fermion 


neutral charged 


neutral 


charged 


neutral 


charged 


y 


O.LO. 2,0.3 


0.1/0.4 


0.L0.4 


0.4 


sei)arati()ii 


200 


500 


250 


# conf. 


400 


600/2,400 


2,400 


800 


3,600 


^ bin 


4 


3/12 


12 


4 


18 


bin size 


100 


200 


200 


source 


point 


wall 



Table 3: 



K 


0.16 


0.145 


size 


8^ 


X 16 


8^ X 16 


12^ X 16 


algorithm 


Metro. 


Cluster 


Metro. 


Cluster 


Metro. 


Ea 


0.809(3) 


0.813(30) 


0.415(2) 


0.419(4) 


0.412(2) 


Zo 


2.78(15) 


2.76(18) 


3.36(13) 


3.33(19) 


3.37(13) 


-E'tt 








0.416(3) 


0.409(4) 


0.410(2) 


Zt^ 


5.94(5) 


5.87(8) 


3.39(17) 


3.34(24) 


3.37(12) 








0.415(1) 


0.414(3) 


0.411(2) 


z, 






6.73(13) 


6.66(23) 


6.75(15) 


V 


0.401(1) 


0.402(1) 










0.303(1) 


0.304(1) 


0.177(1) 


0.177(4) 


0.176(1) 



Table 4: 



y 




0.1 


0.2 


0.3 


En 





0.15(1) 


0.277(7) 


0.384(7) 


E^ 





1.7(2) 


1.76(9) 


1.71(8) 


End 


6 


1.4(3) 


1.3(2) 


1.5(2) 


Ec 


3 


0.149(24) 


0.277(1) 


0.386(16) 




4 


0.152(49) 


0.276(11) 


0.38(2) 




5 


0.168(83) 


0.283(19) 


0.38(7) 



Table 5: 



t 


n 


c 


L 

C 




q3 1 c 

o X Id 


O X ID 


iz X ID 


U 


.yuo(^yj 


n c\QO(ciQ\ /' \/ 1 n— 4\ 
U.Uodl^yoj (^XiU j 


n ^/lQ/'/l^^ /"win— 4\ 
U.U48(^4yj (^XiU j 


i 




0.i6y4(ij 


loo. 6(^3 ) 


Z 


0. 400(^1 j 


4i. 


33.1(^3 j 


o 
O 


z.oll^l ) 


lo.6(D j 


y.3(2 ) 


4 


1.55(1) 




2 6f7l 


5 


1.07(1) 


2.2(1) 


0.88(5) 


6 


7.75(1) 


0.90(1) 


0.33(3) 


7 


6.16(1) 


0.49(6) 


0.15(2) 


8 


5.66(2) 


0.36(8) 


0.12(1) 


9 


6.16(1) 


0.47(9) 


0.17(2) 


10 


7.75(1) 


0.85(1) 


0.35(3) 


11 


1.07(1) 


1.9(1) 


0.87(5) 


12 


1.55(1) 


4.8(3) 


2.8(1) 


13 


2.31(1) 


13.1(5) 


9.3(2) 


14 


3.466(1) 


40.2(6) 


33.1(3) 


15 


5.227(9) 


168.8(1) 


159.9(5) 



Table 6: 



size 


y = 0.1 0.4 


8^ X 16 

En 


0.1895(1) 0.5491(6) 
1.109(1) 1.458(2) 


12^ X 16 E]^ 

Et- 


0.5495(1) 

1.1562(9) 




Figure 1: 
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Figure 4: 
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Figure 7: 
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Figure 8: 




Figure 9: 
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Figure 12: 
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Figure 13: 



